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In this paper we investigate integrable models from the perspective of information theory, exhibit- 
ing various connections. We begin by showing that compressible hydrodynamics for a one-dimesional 
f*"«« ■ isentropic fluid, with an appropriately motivated information theoretic extension, is described by 

a general nonlinear Schrodinger (NLS) equation. Depending on the choice of the enthalpy func- 
tion, one obtains the cubic NLS or other modified NLS equations that have applications in various 
fields. Next, by considering the integrable hierarchy associated with the NLS model, we propose 
^q' higher order information measures which include the Fisher measure as their first member. The 

lowest members of the hiearchy are shown to be included in the expansion of a regularized Kullback- 
Leibler measure while, on the other hand, a suitable combination of the NLS hierarchy leads to a 
Wootters type measure related to a NLS equation with a relativistic dispersion relation. Finally, 
through our approach, we are led to construct an integrable semi-relativistic NLS equation. 
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INTRODUCTION 

Integrable equations are fascinating not just because of their soliton solutions and the connections they make among 
different areas of mathematics, but also because they do describe real physical systems in some limit. An example is 
the cubic nonlinear Schrodinger equation, 

# t +Vx 2; + 2 K 2 |^| 2 V-0 (1) 

which is of relevance in quantum optics, condensed matter physics and other areas. The basic equation |T]) can be 
modified while still preserving integrability, for example by adding to the right-hand-side of ([I]) a term proportional 

06' t0 

R: q= s ^K ( 2 ) 

where p — \ijj\ 2 . Through a change of variables one can actually absorb that extra term and regain the form ([TJ) at the 
expense of redefined parameters However this is possible only if parameter s < 1, whereas if s > 1 one ends up 
with a reaction-diffusion equation . Such Q augmented NLS equations have appeared in plasma physics [1] , where 
they desribe transmission of uni-axial waves in a cold collisionless plasma subject to a transverse magnetic field. 

The reason for using the symbol Q is because such a term, often referred to as a "quantum_potential" , appeared 
first in alternate ways of writing the usual linear Scrodinger equation of quantum mechanics [J [H, Q . Consider the 
one dimensional time-dependent Schrodinger equation (we set the mass m = 1), 

ft 2 

ifupt + y^ xx ~ u ( x )^ = ( 3 ) 
Then substituting into this equation the Madelung representation of the wave function 

i> = (4) 



decomposes it into two real equations, 



St + ^f + U-^-y^ = 0, (5) 
Pt + ( P S X ) X = 0. (6) 
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The first equation may be viewed as a generalisation of the usual Hamilton-Jacobi equation by the term with explicit 
ft dependence, the quantum potential, encoding the quantum aspects of the theory. The second equation is the 
continuity equation expressing the conservation of probability. 

Several attempts have been made to motivate the form of Q and thus obtain a derivation of Schrodinger's equation 
from classical dynamics Q. Here we adopt an information theoretic perspective similar to that used in statistical 
mechanics and which is usually refered to as the "maximum entropy method" Q. The idea is that if one has a 
system that has to be described probabilistically then, lacking any information of the detailed microscopic dynamics, 
one should choose the probability distribution with minimum bias. This is achieved by maximising an appropriate 
measure of uncertainty (entropy), such as the Gibbs-Shannon measure used in classical statistical mechanics. 

In order to proceed with an information theoretic interpretation of (|5j [6]) , it is useful to approach those equations 
through a variational principle [9j : one minimises the action 



$ = 



s t + l(s x ) 2 + u 



dxdt 



-If 



(7) 



with respect to the field variables p and S. The positive quantity 



dxdt p 



Px 



(8) 



resembles the "Fisher information" measure used in statistics [id 111], Since a broader probability distribution p(x) 
represents a greater uncertainty in x, so If may be thought of as an inverse uncertainty measure. 

The equations (JTJ [SJ) were used in Ref.Q to interpret Schrodinger's equation as follows: First one notes that without 
the term Ip , varying Eq. ([7]) gives rise to the Hamilton-Jacobi equation describing a classical ensemble. The probability, 
p(x) appears in this context because one supposes that there is uncertainty in our knowledge of the initial position 
of the particle. One then adopts the principle of maximum uncertainty [8y to constrain the probability distribution 
p(x) characterising the ensemble: we would like to be as unbiased as possible in its choice, consistent with our lack of 
information. That constraint is implemented in |(7J| by minimising Ip when varying the classical action: ft 2 / 8 is the 
Lagrange multiplier. 

It remains to explain why Ip is chosen as the information measure in the above quantum mechanical context as 
opposed to say the Gibbs-Shannon measure. In information theory and statistical mechanics the Gibbs-Shannon 
measure is the simplest possibility that satisfies certain axioms that are deemed necessary in those contexts Q. 
Similarly one can derive the Fisher measure as the relevant quantity that satisfies axioms relevant for classical ensemble 
dynamics and hence appropriate for use in deriving Schrodinger's equation [12j 

In this paper we would like to apply the above information theoretic reasoning to motivate the NLS ([1]) and its 
various extensions. In the next section we first review the derivation of the action for a classical compressible fluid 
in one dimension. Then in Sect. 3 we use information theoretic arguments to modify the action and so arrive at a 
general nonlinear Schrodinger equation. In Sect. 4 we employ an expansion of the enthalpy function to obtain specific 
examples of the nonlinear Schrodinger equation. In Section 5 we consider the hierarchy associated with the NLS 
equation and use that to define a hierarchy of higher-derivative information measures. Then in Section 6 we relate the 
information hiearchy to other information measures in the literature and in Section 7 we use the information measures 
to construct NLS equations with relativistic dispersion relations. Our conclusion is in Section 8 while in the appendix 
we discuss the relationship between the Madelung representation of quantum mechanics and the complexified Burgers 
equation. 



COMPRESSIBLE FLUID IN ONE DIMENSION 



In this paper we focus on a specific physical model, hydrodynamics, to illustrate our approach, though we believe 
that much of it can be generalised to other contexts. 

The Euler and continuity equations for a one-dimensional compressible fluid are 



-P x = 0. 
P 



(9) 



Pt + (pv) x = 



(10) 
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where the hydrodynamical variables corresponding to the density of fluid, velocity of fluid and pressure have been 
denoted by p[x, t), v(x, i) and P(x,t) respectively. In addition one has the thermodynamic equation [l3[ 

_ = _+„§, = () (11) 

expressing the conservation of entropy, §, if one assumes the absence of the heat exchange between parts of the 
medium. To complete the description of the dynamics one also needs an equation of state 

P = P(lp) (12) 

whose concrete form depends on the properties of the fluid. For example, an ideal gas has 

P = e § / Cl V , (13) 

where 7 = cp/cy (the Poisson adiabate) is the ratio of specific heat capacities at constant pressure and volume 
respectively. 

We will consider an isentropic fluid which has a spacetime constant § so that (fTTj) is automatically satisfied. For 
such barotropic processes the pressure becomes a function of density only, 

P = P{p) . (14) 

It is convenient to introduce the enthalpy function defined through the relation 

which implies, for p x =/= 0, 

E\p) = -P'(p) (16) 



or 



P 



rp dP 

E{p)= I (17) 



Po 



p 



Then the system of equations (|9"j) . (fTU|) becomes 

v t + vv x + (E(p)) x = 0, (18) 

Pt + {pv) x =0. (19) 
This system may be written in Lagrangian form by first introducing the velocity potential 

v(x,t) = S x (x,t) (20) 
then integrating the first equation once and introducing the enthalpy potential 

•m - ^ (2D 

dp 



to get 



The action for this system is 



Pt + {pS x ) x = 0. (23) 



A = I [pS t + + v (p) ) dxdt ( 24 ) 



and equations (|22|) . (f23|) appear by varying this functional with respect to p and S. The resemblance of (|24)) to the 
classical part of (O will be the starting point for the extension in the next section. 

We note, for later use below, that when the enthalpy vanishes, E — 0, the fluid equations (|22| . ([23|) are invariant 
under a scaling, p — > ap. That is, with E = 0, the magnitude of the density does not matter, only its variation. 
When E is not zero, the equations to be derived later become generalised nonlinear Schrodinger equations, with a 
sensitivity to the magnitude of p. 



4 



INFORMATION-THEORETICAL EXTENSION OF COMPRESSIBLE FLUID DYNAMICS 



The action ([24]) gives the classical equations of motion for the fluid. Since the density p informs us about the 
likelihood of finding the microscopic fluid elements at a certain region of spacetime, it plays a role analogous to 
the probability density in quantum mechanics. Thus from the density we may form an information measure I that 
quantifies our knowledge of the microstates and we may demand, as in the previous section, that the equations of 
motion follow from (|24p but constrained such that our uncertainty (information) is maximised (minimised) . This will 
lead to modified hydrodynamics equations that depend on the form of information measure chosen in the procedure. 

Now, the density p(x) is positive definite and if it is uniform it tells us that the underlying particles of the fluid 
could be anywhere: we have no information (maximum uncertainty). If the density is peaked somewhere, we know 
that a fluid particle is more likely to be there, that is we have gained information. Thus we require that our scalar 
information functional I[p] have the property that it is positive definite and / — » as p — > a constant. 

We prefer local equations, and so we may write I as an integral over a density function J(p), 

dxdt pj(p) (25) 

Next we assume the density to be slowly varying and so do a derivative expansion, 

J(p) = J {p) + p'Ji(p) + p".h\(p) + {p'fJ22{p) + higher derivative terms , (26) 

where Jo(p), Ji(p), ... do not contain any derivatives. We assume that when (|26[) is used in (|25p the integrals are 
convergent term by term. 

We also impose the strong condition that the information measure, /, does not break the invariance of the E = 
equations of motion (f22|) and (|23|) under the scaling of p — > ap. That is, although E will generally break that 
invariance, we demand that the terms in the modified equations of motion that come from / do not do so: the 
information measure is chosen to be neutral to the magnitude of p but measures only local variations. So here we see 
the first difference between the contributions of our / and E : I is unbiased towards the size of p. 

In order to achieve our goal, we need to demand that J{p) in (|25[) is scale invariant (note we already factored out 
a p in the integral form of I). In that way the equation of motion terms that come from varying / will be scale 
invariant. This is satisfied if (|26|) has the form 

J{p) = a + Oi x (p'/p) + oai x (p"/p) + a 22 x (p'/p) 2 + , (27) 

where a*,, aui, • ■ for k, I = 1, 2, ... are constants. Using this in the integral that defines /, Eq. (|2"5"l) . and dropping constants 
and total derivatives, only the 022 term survives to leading order and it gives precisely the Fisher information measure! 

Recall that we still need to demand positivity of our information measure: that fixes the lagrange multiplier to be 
positive if we are minimizing the information. Fortunately, the Fisher measure already satisfies the other required 
property, that it vanishes as p — *■ a constant. 

Note that in the fluid problem we have in general p — > po as \x\ — > 00, so that 



(p(x,t)-po)dx = l. (28) 

With the choice of the Fisher information measure 

If = J ^J^ dx = 4 / (VpUVphdx (29) 

as motivated above, the extension of the boundary condition l|28j) from the usual case in quantum mechanics does not 
modify the convergence properties of I. 
Thus we have variational functional 

a +y If (30) 

where A is the Lagrange multiplier and the equations of motion that follow are 

(S x f dV{p) A 2 (^p) xx _ 

St + — + — -Y^p^-°> (31) 
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Pt + {pS x ) x = 0. (32) 

These equations may be combined into one complex equation through the inverse Madelung transformation to give 
the following general Nonlinear Schrodinger equation, 

iA^ + y^x--B(|V| 2 )V = 0. (33) 

In summary, although both the enthalpy function, E, and the information functional, /, will contribute p dependent 
terms to the equations of motion, their structure and origin is in general different. Using the lowest order information 
measure with the properties described above we get a generalized NLS equation (|33[) that depends on the form of E. 
Using more generalised information measures will give further extensions of the NLS equations. 

We remark also that unlike the quantum mechanics case [T^, the deduction of the Fisher measure above did not 
use the separability condition: rather here we assumed a convergent derivative expansion of the information density. 
See also [14J for similar arguments used in the relativistic case. 



INTEGRABLE CASES 



In this section we study the simplest form of the function E in (|33)) that will give rise to integrable systems. If E 
as function of p = \ip\ 2 is analytic then 

£(|V| 2 ) = Eo + EM 2 + EM 4 + - + £„M 2 " + - (34) 
This equation implies that the pressure according to (|15p is also an analytic function of the form 

1 o n 

P(p) =E Q + -E lP 2 + -E 2 p 3 + ...—-E n p n+1 + ... (35) 

2 3 n+1 

At the lowest order of nonlinearity we get the Nonlinear Schrodinger Equation (NLS) with cubic nonlinearity 

zA^ t + y il> xx - {E + Bk | Vf> = (36) 

This model is integrable for both signs of E\. For E\ > it is defocusing (e.g. repulsive Bose gas) and nontrivial 
soliton solutions exist only with nontrivial boundary conditions, so in this case Eg Q. For E\ < we have the 
focusing case (e.g. attractive Bose gas) for which soliton solutions exist for vanishing boundary conditions; so in this 
case we can put Eq = 0. Both cases have applications in nonlinear optics describing pulse propagating in nonlinear 
media. Consider the second case again: by rescaling space and time variables t' = t/X, x' = \/2x/\ and the coupling 
constant E\ — —2k 2 we may rewrite it in the form (we now skip all upperscripts) 

^t + ^ + 2 K 2 |^|V = 0. (37) 

The "semiclassical" or dispersionless limit of this equation was studied in [l5| in relation to shock wave propagation 
in nonlinear optics. The wave form of this semiclassical limit is a NLS equation perturbed by a quantum potential 



Ivl 

Then we can conclude that inclusion of information characteristics in the form of the Fisher measure, produces NLS 
(157)1 from dispersionless NLS ([3"5]) and corresponding solitons of the first one from the shock waves of the second one 



INTEGRABLE NLS HIERARCHY AND HIGHER-DERIVATIVE INFORMATION MEASURES 



It is well-known that one can construct a hierarchy of higher-order differential equations that are related to the 
cubic NLS (|3~T[) and its complex conjugate, which are still integrable 17j . 



i a 3 



= K- 



AT ft 



(39) 
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where tjy, N = 1,2, 3, ... is an infinite time hierarchy. Here TZ is the matrix integro-differential operator - the recursion 
operator of the NLS hierarchy - 

/ d x + 2«V f x i> -2k 2 iIjJ x iP \ 
K = ia 3 [ . (40) 

\ ~2k 2 ?PJ x iP d x + 2K 2 ^J x i;J 

For the first few members of the hierarchy N = 1,2,3,4 this gives 

= V>* , (41) 



tyfa +^+2k 2 |^|V = 0, (42) 



^t. +^ + 6k 2 |^| 2 ^ = 0, (43) 



#t 4 = Vw* + 2k 2 (2|^| 2 V' + 4|V|Vx* + V^ 2 + 3^V 2 ) + 6k 4 |V|V (44) 
In the linear approximation, when k = 0, the recursion operator is just the momentum operator 

Wo=iff 3 ^ (45) 

and the NLS hierarchy (|39p becomes the linear Schrodinger hierarchy 

i^ tn = i n d^ . (46) 

Written in the Madelung representation it produces the complex Burgers hierarchy (see Appendix for n — 2 case and 
[H | for arbitrary n). 

Let us look more explicitly at the fourth order flow for which the Hamiltonian is 

H = J [i> xx ^ xx - 8k 2 ^ x ^ - K 2 {i, 2 x i> 2 + ^ 2 X ) + 2k 4 \^\ 6 ] dx . (47) 
In the Madelung representation 

4, = ^pe lS = e R+lS (48) 

this becomes 

H f rPL _ pxxpI ji_ (Ag) 

J Up 4 P 2 16p 3 + 1 ' 



P St x + pS x + 2p x S x S xx - 2k' ^-pi + 3p'S z x ) + 2K*p A ]dx (50) 

In fact for configurations with S = const we have only contributions from the first three terms which can be combined 
into 

H = I (\fp)xx{\fp)xxdx . (51) 



This may be considered as a higher order analog of the Fisher information measure If (j29|) . 

Generalizing, the above linearized Schrodinger hierarchy suggests, after the substitution S = const, the even order 
information measure hierarchy 



h = J {y/p)x{yfp)xdx , 



(52) 
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(y/p)xx(y/p)xxdx , (53) 



hn = J (y/p)x..x(y/p)x...xdx, (54) 

Here all odd members vanish because their integrands are total derivatives. 

We will use the above information hierarchy in the next section to construct relativistic NLS equations and exhibit 
links between different information measures known in the literature. 

Before leaving this section, we compute the above information measures for the one soliton solution of the NLS 
equation. The measure is 

P"( x ) = n 72 ( 55 ) 

2 cosh vx 

which satisfies 

/>oo 

p v (x)dx = 1 (56) 



and is characterized by a real parameter v so that 

lim p v {x) = S(x) . (57) 

v — >oo 

Then information measures (]52H . (|53p.... for this distribution apart from numerical constants are simply 

h=v 2 , h = v\ hu = v 2n ,.... (58) 

RELATION TO KULLBACK-LIEBLER AND WOOTERS MEASURES 

The Gibbs-Shannon entropy 

I GS = - J p{x) In p(x) dx (59) 



may be genelaised to the Kullback-Leibler information 11 1 



I KL (p,r)= - [ p(x)ln^-dx, (60) 
J r(x) 

where r(x) is a reference probability distribution. If one chooses the reference distribution to be the same as p(x) but 
with infinitesimally shifted arguments, that is r(x) — p(x + Ax), then to lowest order, 



I KL (p(x),p(x + A(x)) = J^Ll F (p( x )) + 0(Ax) 3 , (61) 

that is, the Fisher measure is recovered as the lowest order term in the expansion. 

One may further generalise the Kullback-Liebler information by introducing a parameter < r] < 1, as used for 
example in [l9| . 

f p ( x ^ n (i -rf)p(x) + r]p(x + i]L) dX ' ^ 
This form is nonsingular even if the density vanishes at any point. For L << 1 we have the expansion 

06 „5 r „3 
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rf 
24 



Pxx J , / V 

dx+ (— 

P 3 



-) 



0(L 5 ) 



where a number of surface terms have been dropped after integration by parts. Let us look at the symmetrised 
measure 



M(+L) + M(-L) = iV / &dx- 

J P 



—L 



12 J P l 3 9 A' J p\ 



0(L 6 ) 



where as before, the lowest order term, proportional to L 2 , is the Fisher measure. By choosing the parameter -q to 
satisfy 



2 4 3 „ 

T) — —T) + — = , 



(63) 



or r\ = (2 ± *\/5/8)/3 we can rewrite the next, 0(L A ) term, as the higher-derivative information measure given by 

Thus the first two members of the information hierachy l2 n we proposed in Sect. (5) are contained in the Kullback- 
Leiblcr information. 



RELATIVISTIC NLS EQUATIONS 

Now let us consider other ways of combining and using the information measures. Take a Hamiltonian of the form 

H = C 2 h + C 4 I 4 + •■• + C 2n hn + •■■ (64) 

where the constant coefficients Cj depend on the context. For example, for low momenta one may expand the relativistic 
dispersion relation E = y m 2 c 4 + p 2 c 2 to get 



E — mc + - — 



P 



2m 8m 3 c 2 

This may be used to construct a "semi-relativistic" Schrodinger equation as a formal power series 

h 4 d 4 



1 dV mC [ 2m 2 c 2 dx 2 8m 4 c 4 5a; 4 
= hip 



(65) 

(66) 
(67) 



where apart from a constant, the average of h for real ip is precisely for a particular choice of coefficients. In fact 
one may proceed further: by replacing the derivative operator d/dx with the recursion operator 1Z one obtains an 
integrable relativistic nonlinear Schrodinger equation 



ia 3 



mc 



,n 2 



(68) 



We note that relativistic versions of the Schrodinger equation have been considered in different contexts, for example 
to study relativistic quarks in nuclei [2fj] and gravitational collapse of a boson star 21 1. A nonlinear version has 
appeared as the semi-relativistic Hartree-Fock equation [22] |. But none of those models is known to be integrable. 
By contrast the model (|68p . where the square root is considered as a formal power series (matrix pseudo-differential 
operator), is an integrable nonlinear Schrodinger equation with relativistic dispersion: 



mc 



1 



d 2 



m 2 c 2 dx 2 



(69) 



9 



where the nonlinearity expanded in 1/c is the infinite sum 

FW = ^[-2k 2 MV] 



8m 3 c 2 

(70) 

What is amazing is that if we expand also the dispersion part in 1/c 2 , then at every order of 1/c 2 we get an 
integrable system. It means that we have integrable relativistic corrections to the NLS equation at any order. 

Another way of constructing a relativistic model that includes higher-derivative information measures is to use 
rapidity variables for the relativistic dispersion relation, 

E = mc 2 coshx, p = mcsvahx , (71) 

This gives the relativistic model with Hamiltonian [23| 

f d 

H = mc 2 ip cosh(L-) tfjdx (72) 

Expanding the cosh we again have a member of the information hierarchy and a relativistic NLS equation. 
Finally, the above free Hamiltonian may be represented as a finite difference operator 

H = mc 2 /2 / ^(e L £ + e- L &)ipdx = / (^(x)ip(x + L) + ^{x^x - L))dx . (73) 



The dispersive part of this hierarchy for S — const gives a Wootters type [24j measure 



!w = / (y/plxjy/p(x + L) + yfpjx)yj p{x - L))dx . (74) 



SUMMARY 

We have shown how information theory arguments can be used to motivate the general nonlinear Schrodinger 
equation in the context of hydrodynamics. This then led us to study different information measures. 

We noted that the integrable hierarchy of linear and nonlinear Schrodinger equations, in their Madelung form, 
naturally suggest a hierarchy of information measures of which the Fisher measure represents the first member. The 
lowest members of the information hiearchy were shown to be included in the expansion of a regularized Kullback- 
Leibler measure. 

We also showed that how to contruct integrable semi-rclativistic nonlinear Schrodinger equations using various 



combinations of the information measures. These classes of equations, which are distinct from those obtained in [14 1 
and references therein, might be useful in analyzing relativistic corrections to solitons, Bose-Einstein condensates or 
other condensed matter systems with effective equations of relativistic form. 

Acknowledgment This work was supported in part by the National University of Singapore and the Izmir Institute 
of Technology, Turkey. 

Appendix A: Madelung Fluid and Complex Burgers Equation 

Here we show that the Madelung representation used in quantum mechanics may be viewed as a complexified 



version of the Cole-Hopf transformation 25j, [26| which relates the linear heat equation with the nonlinear Burgers 
equation [27j |. 

Begin with the one dimensional time-dependent Schrodinger equation, 

h 2 

ihipt + 7T~ Vfe* - u ( x )^ = ( 75 ) 
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and define complex velocity 



we have 



u(x,t) = u c + iu q = - — (ln^)aj (76) 



u c = ^S x , u q = -^-(lnp) x {,,) 



This show that the real part of the complex velocity is the classical velocity, while the imaginary part can be associated 
with the " quantum velocity" 28j ■ Considering (f76|) as complex analog of the Cole-Hopf relation we can see that the 
Schrodingcr equation (|T5[) is equivalent to the complex Burgers equation 

ih 1 

u t + uu x = —u xx U (78) 

Im m 

Splitting this equation into real and imaginary parts, for the former we have a hydrodynamic equation for the classical 
velocity 



1 ,„ n 



rn 



u ct +u c u cx = - — {U+-u qx -—u q ) x (79) 
while for the imaginary part we obtain the transport equation for the quantum velocity 

u qt + (u c u q ) x = -^u cx (80) 
After one space integration the last equaton becomes the continuity equation 

Pt + {pu c ) x =0. (81) 

Equations ((79|) .(I81 |) describe a Madelung fluid. Thus the above consideration show that a Madclung fluid can be 
described by one nonlinear complex Burgers equation (|75j) . In equation (|79")l the potential part defines a quantum 
potential written through the Riccati equation 



h m 

U q = ^u qx - -u q 2 . (82) 



Then expression (|77|) for the quantum velocity 

?i Px 



gives a linearization of the Riccati equation 



(83) 
Am p 



2m 



with potential 



ft 2 (y/p)xx _ ft 2 M XX 

2m yfp 2m 



TT _ i" \\ F)xx _ irixx , at .\ 

U 1- ~ ~ -~~~UA~- 



Finally we note that equation (184|) has simple geometrical interpretation. If 

ds 2 = dt 2 + \i)\ 2 dx 2 (86) 
is the metric of a surface with geodesic coordinates (t,x), then the Gaussian curvature becomes the quantum potential 



( VP)xx _ 2m, 
VP " Ti 2 



K=- y -^± = ^lu q (87) 
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This shows that the quantum potential can be interpreted as the Gaussian curvature 

U q = ^-K. (88) 
2m 

In particularly it implies that the classical limit corresponding to a vanishing quantum potential locally is a plane. 
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